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We study a correlation function for the one-dimensional isotropic XY model {XXQ 
model), which is called the Emptiness Formation Probability (EFP). It is the probability 
of the formation of a ferromagnetic string in the anti-ferromagnetic ground state. Using 
the expression of the EFP as a Toeplitz determinant, we discuss its asymptotic behaviors. 
We also compare the analytical results with numerical calculations as the density-matrix 
renormalization group and the quantum Monte-Carlo method. 
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§1. Introduction 

Recently there have been a lot of progresses for the calculation of correlation functions for the 
one-dimensional XXZ model, ^"^^ 

TV N 

i=i i=i 

which is exactly solved by the Bethe ansatz method. Among various correlation functions, the so 
called emptiness formation probability (EFP) P{n) 

Pin) = {fl{S] + h), (1.2) 

is the most fundamental one.'^^ 

Concerning the asymptotic behaviours of P{n) for |A| < 1, the following properties have been 
discussed: 

(1) P(n) decays like a Gaussian at zero temperature.^^ 

(2) P{n) decays exponentially at finite temperature. ^-^ 

In this paper, we shall confirm these properties for the isotropic XY model (A = 0)^'^'^^-' in detail. 
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Actually, based on the representation of P(n) as a Predholm determinant, (c.f. eq. (2.13)), the 
property (1) has been established for generic |A| < 1 to some extent . 

For the special case, A = 0, we will show in this paper that P{n) can be alternatively repre- 
sented by a Toeplitz determinant. The representation by Toeplitz determinant is very powerful to 
give us more information on the properties of P(n). Actually there are theorems known for the 
asymptotics of the Toeplitz determinants, which directly establish properties (1) and (2). They 
also give us the complete forms of the asymptotic behaviours including the prefactor. 

On the other hand, we can calculate P{n) for small n by other numerical methods. In this 
paper we employ the density-matrix renormalization group (DMRG) method for zero temperature 
case and the quantum Monte-Carlo (QMC) method for finite temperature case. The numerical 
data calculated by the DMRG and the QMC agree satisfactorily with the analytical ones. We 
believe the method and the results presented in this paper will afford a good foundation when we 
investigate the EFP for other models. 

The plan of the paper is as follows. In §2, we derive the representation of P{n) by Toeplitz 
determinant. The asymptotic form of P(n) is obtained from Widom's theorem. In §3, we demon- 
strate a numerical simulation of P{n) by DMRG and compare with the analytical calculation. In 
§4, the EFP at the finite temperature is studied. This time the asymptotic function is obtained 
from Szego's theorem. It is compared with the numerical calculations by QMC in §5. The last 
section is devoted to summary and discussions. 

§2. EFP at Zero Temperature 

The Hamiltonian of the isotropic XY model is defined by the special case of (1.1) with A = 0, 



I.e., 



N N 



where are the spin 1/2 operators 



[S'^,5f] — iSkieap^Sl, 



(2.1) 



(2.2) 



and are represented by the Pauli matrices as 



^3 ~ 2 ^' ' 



0"o 



"o 


f— ^ 


V 


"O -i" 




"l " 




'-1 = 


1 





_i _ 


-1_ 



In the Hamiltonian (2.1), we assume the periodic boundary condition (PBC), 

OCX QOi 



(2.3) 



(2.4) 
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We also introduce 5^ = ± iS^ for later convenience. The EFP at zero temperature is defined 

by 

n ^ 

P(n) = {GS\ lliS] + -)\GS), (2.5) 

i=i 

which represents the probability of the formation of a ferromagnetic string with a length n in the 
anti- ferromagnetic ground state \GS). 

The Hamiltonian (2.1) can be diagonalized by means of the Jordan- Wigner transformation/'^) 

fe-i 

fe-i 

S^ = ll('^-^cj)ck, (2.6) 

where c|.,Cjk are the fermion creation and annihilation operators, respectively. In fact, the Hamil- 
tonian (2.1) is transformed to a spinless fermion model as 

N-l N 

2 



3=1 j=l 

+ ^ 1 11(1 - 24cfe)| {c\cN + 4ci). (2.7) 

We can diagonalizc the Hamiltonian (2.7) by means of the Fourier transformation. In the thermo- 
dynamic limit (A^ — > oo), the ground state energy is given by 

1 r''^ 

e = -h+— {2h- Jcosq)dq, (2.8) 

where the Fermi point is determined by 

fep = cos-^(2/i/J). (2.9) 

Note that kp goes to zero ash^hc = J/2. The EFP (2.5) can be calculated by use of the Fermion 
operators as 

n 

P{n) = iUcjc]) 
i=i 

= det[(c/cj„)]]^„=i 



-I n 

det 



27r 



(2.10) 

J l,m=l 

where we have used Wick's theorem. A kind of the determinant like eq. (2.10) is referred to as a 
Toeplitz determinant. 
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It is known that the transverse correlation function {S^Sq) for the XY model is also written 
as a slightly different Toeplitz determinant. -"^^^ It decays algebraically as (r — oo), 

{S^^SS) - l/V^. (2.11) 

Similarly, the longitudinal correlation function decays algebraically 

{{S^^S^)) ^ {S^^S^) - {S^){S^) ~ l/r'. (2.12) 

For further results on these correlation functions, see the refs. 10,12-20. 

As we shall show below, the asymptotic behavior of P{n) is remarkably different from the two 
point correlation functions above. 

We remark that there are two other representations for P{n). 
(1) Predholm determinant representation:^' -"^^^ 

P(n) = det [I - Kn] , (2.13) 

where is a linear integral operator acting on functions on interval [— fcp, k-p] with a kernel 

I n— 1 



Or equivalently. 



27r 

e 2 ' sm ^ (p - g) 
27r sin^(p — g) 



(2.14) 



^(") = E —r / dgi • • • / dg,- det [Kn{quqm)]im=i 

j=0 -/ -/i^F J -k-p 

/'27r— /uF /•27r— /cp 

= -r / dgi • • • / dg„ det [i^n(g/, 9m)]rm=i • (2-15) 



^27r— /uF /•27r— /cp 

dgi--- / 

' fcp J fcp 

For the derivation of the equalities above, we refer to rcf. 24 
Note that the kernel Kn{p,q) (2.14) can be replaced by 

sin§(p-g) 
27r sin ^ (p — g) ' 

due to a property of the determinant.^^) 
(2) Multiple integral representation:^' 



(2.16) 



s t «(«-!) n(„+i) f dAi [ dA„ T-r 
P{n) = (-1) 2 2 2 / ... / II 



sinh2(Aa - Afe) 
sinh(Aa - Afe - §i) 



^^f^ si^h^"HA,■-f^)sinh"-^^(A, + f^) 

""-Li cosh"2A, ■ ^^-^^^ 

j=i J 
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Here the integration range R is defined by 

r-A 



dA 

R 



dA for A;f > f 

A 

A 



/ + / + / dA for fcp < f 

J —oo J —oo J A / 

where the endpoint A is related to kp through 



(2.18) 



F 



e-2A = -i^,^ 1. (2.19) 

gifcp _^ i ^ J 

The integral representation (2.17) was first obtained in ref. 25, in the case of A = oo. Recently it 
was confirmed and generalized to finite A case in ref. 6. 

We can show that eqs. (2.10), (2.13) and (2.17), are ah equivalent.^^'^"^'^^) For the readers' 
convenience, we give some brief proofs in Appendix A and Appendix B. 

Note that the Fredholm determinant representation (1) and the multiple integral representation 
(2) do exist for the general XXZ model (1.1),^'^"^'^^) while the Toeplitz determinant (2.10) is very 
special to the XY model. However, in the following, we discuss the asymptotic properties of P{n) 
for the isotropic XY model (2.1) based on the representation (2.10). 

Let us first assume that h takes some moderate value, \h\ < he- Then from Widom's theorem^'''-* 
for the Toeplitz determinant (2.10), we immediately obtain 

In P(n) = In cos y - ^ ln(n sin y ) + ^ In 2 + 3C' (- 1) + o(l), (2.20) 

i.e., 

1 2 

P{n) ~ 2iie3^'(-i) (sin^^ ' n-V^ (cos , (n ^ oo), (2.21) 

where ("(x) is Riemann's zeta function. We note that (^'(—1) = —0.165421 • • • and 

2i^e^^'(-^) = 0.645002 • • • . (2.22) 
In terms of the external magnetic field h, the formula (2.21) can be expressed as 



1 "2 



Pin) . 2^e3C'(-i) ("-f/') n-i (^^) \ (n ^ oo). (2.23) 

Thus P{n) for the isotropic XY model (2.1) decays like a Gaussian. The speed of Gaussian decay 
as well as the prcfactor depends on the magnetic field h. The fact that P{n) shows Gaussian 
decay has been established in ref. 5, 23 and 26. Especially wc remark that the dependence on the 
magnetic field h in the part of Gaussian decay was already obtained in ref. 5. Their method is 
based on the solution of the Riemann-Hilbert problem associated with the Fredholm determinant 
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(2.13). However, the prefactor and the part n~^/^ have not been obtained from the Riemann-Hilbert 
approach yet. 

We observe that the prefactor of (2.20) diverges as h gets closer to the critical value, i.e., 
h ^ hc = J/2. There the asymptotic formula (2.23) may not be adequate. For this region, as /cp 
is very small, we can expand the determinant (2.10) with respect to A;f. Actually we can easily 
obtain the expansion. 



(2.24) 



which may be useful to study the critical behaviors of P{n) near h = h^- 

In Fig. 1, we plot P{n) and the corresponding asymptotic function (2.21) for several magneti- 
zations m = (<S'|). Here note that m is connected to as 



1 A:f 
m = . 

2 vr 

The readers will observe P{n) is well fitted by the asymptotic function (2.21). 



(2.25) 
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Fig. 1. P(n) at zero temperature. 

We would like to remark that the Toeplitz determinant (2.10) appears in relation to the level 
spacing distribution of the random matrix theory. More precisely, the function ^2(0) for Dyson's 
circular ensemble of n x n unitary matrices, is identical to eq. (2.10) with the correspondence 
a = fcp. Here ^2(0:) means the probability that an interval (of the unit circle) of length 2a contains 
no eigenvalues. The function A2{a) has been studied profoundly in the context of the random 
matrix theory. Then we can translate all the results obtained for ^2(0) to P{n). For example, 
it is recently shown that ^2(0;) satisfies an integrable non-linear differential equation, which is 
equivalent to Painleve VI equation.^^'^^^ Hence our P{n) also satisfies the Painleve VI equation.^^^ 
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§3. Numerical Simulation by DMRG Method 

In this section, we perform numerical simulation so as to confirm the validity of the asymp- 
totic form (2.21). We employed the density-matrix renormalization group (DMRG) here.2^'29)The 
method has an advantage in that it allows us to treat large system sizes. Our algorithm here is 
standard, and its detail would be found in literatures; we refer the readers to ref. 30. Hence, 
below, we will outline some technical points that are relevant to our simulation precision: We im- 
plemented the infinite-system method, which is adequate to study the ground-state properties in 
thermodynamic limit. We have repeated hundreds of renornializations. At each renormalization, 
we remained, at most, two- hundred relevant states (bases) for a block; namely, in conventional 
terminology, we set m = 200. The density-matrix eigenvalue {wa} of each remained base indicates 
its significance (weight). We found Wa > 10"^*^. That is, weights of discarded states are no more 
than Wa ~ 10~^°. 

However, there are some subtleties in the DMRG calculation when S = 1/2 as in our present 
case. The number of spins consisting each block increases by one after another through each 
renormalization. The problem is that the structure of the Hilbert space changes significantly 
according to the situations whether block contains even number of spins or odd number of spins. 
(Note that for S = I, such difficulty does not arise.) Therefore simulation data alternate in turn 
through renormalizations, even though we repeat hundreds of renormalizations. In this respect, the 
translational symmetry is broken intrinsically in the scheme. In other words, edge effect remains 
through renormalizations. In fact, as for P{2), for instance, the digit of order 10~^ alternates. 
We coped with the difficulty by averaging the data over two successive renormalizations. Then 
we could achieve the precision of 10~^ for P{2). The precision of data would be improved further 
if we Tisc a certain careful extrapolation. However, in the following, we shall show that our data 
reproduces the analytical results rather satisfactorily. 

In Fig. 2, we plot the data for P{n) obtained by the DMRG simulation. The dotted lines 
represent the asymptotic functions (2.21). For the positive magnetization m = 0,0.2,0.4, we also 
list the explicit numerical data in Table I-III. We find our numerical data by DMRG agree well 
with the analytical ones in §2. 

We see that for very small P{n) < 10~^^, simulation data start to deviate widely from the 
analytical ones. This may be simply due to the numerical round-off error; note that (double 
precision) real number is stored as sixteen digits in computer and the data of very small P{n) are 
not reliable in principle. For large /i(~ 0.4), P{n) decays slowly and our data yield results reliable 
even for very large distances. Actually, from Table 3, we see that the deviation is less than few 
percents even for n = 10. 
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Fig. 2. Numerical calculation of P{n) by DMRG. 










Table I. P{n) for (S^) = 








n 


DMRG 








exact values 


1 


0.500 








0.5 


2 


1.486 X 10" 


-1 






1.487 X 10-1 


3 


2.367 X 10" 


-2 






2.368 X 10-2 


4 


1.940 X 10" 


-3 






1.945 X 10-3 


5 


8.096 X 10" 


-5 






8.126 X 10-^ 


6 


1.699 X 10" 


-6 






1.715 X 10-*^ 


7 

8 


1.795 X 10" 
9.559 X 10" 


-8 
11 






1.823 X 10-^ 
9.740 X 10-11 


9 


4.578 X 10" 


13 






2.612 X 10-13 


10 


2,442 X 10" 


14 






3.514 X 10-1^ 
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Table II. P(n) for = 0.2 



9 



n 


DMRG 


exact values 


1 


0.699 


0.7 


o 
z 


4.ziy X iu 


4.ZO/ X iU 


3 


2.123 X 10-2 


2.140 X 10-2 


4 


8.717 X 10-2 


8.840 X 10-2 


5 


2.890 X 10-2 


2.950 X 10-2 


6 


7.673 X 10-3 


7.907 X 10-3 


7 


1.628 X 10-3 


1.695 X 10-3 


8 


2.750 X 10-^ 


2.902 X 10-4 


9 


3.713 X 10-5 


3.960 X 10-5 


10 


4.004 X 10-6 


4.304 X 10-6 


Table III. P(n) for {S]) = 0.4 


n 


DMRG 


exact values 


1 


0.899 


0.9 


2 


7.984 X 10-1 


8.003 X 10-1 


3 


6.993 X 10-1 


7.019 X 10-1 


4 


6.028 X 10-1 


6.061 X 10-1 


5 


5.103 X 10-1 


5.146 X 10-1 


6 


4.238 X 10-1 


4.290 X 10-1 


7 


3.455 X 10-1 


3.508 X 10-1 


8 


2.758 X 10-1 


2.811 X 10-1 


9 


2.161 X 10-1 


2.206 X 10-1 


10 


1.655 X 10-1 


1.695 X 10-1 
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Prom the numerical data of P{n) for small n, we can estimate a plausible asymptotic form, 

P(n) ~ const X a"'"' n". (3.1) 

In fact, by fitting successive three data for P{n) in the form (3.1), we get the approximate values 
for const, a and a. We show an example in Table IV, where each parameter is obtained from 
P{k),P{k + 1) and P{k + 2). The exact parameters in the asymptotic function (2.21) are 



Table IV. Estimation of the fitting parameters by DMRG: = 



k 


a 


const 


a 


1 


1.415 


0.707 


-0.247 


2 


1.415 


0.709 


-0.251 


3 


1.414 


0.719 


-0.271 


4 


1.416 


0.673 


-0.204 


5 


1.415 


0.723 


-0.265 



a = 1.414, const = 0.703, a = -0.25. 



(3.2) 



Thus we could perfectly estimate the speed of Gaussian decay, a. We also got reasonable approxi- 
mate values for the pref actor and the exponent a. Prom these results, we can insist that the DMRG 
method provides a powerful means to study the EPF at zero temperature for the ID spin systems. 



§4. EFP at Finite Temperature 

At finite temperature (T > 0), the EPF is defined by 



Pin) = iUiS! + ^))t ^ 

3=1 



Tr(e-«/T) 



(4.1) 



Por the isotropic XY model, we can represent eq. (4.1) in terms of a Toeplitz determinant by use 
of the finite temperature version of Wick's theorem (Bloch-de Dominicis theorem) as follows. 



P(n) = det {cicI^)t 



n 

l,m=l 



det 



{l—m)q 



p{q)dq 



(4.2) 



',,m=l 
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where 

1 1 



1 + e-^(«)/^ 1 + e<i)/'^ ' 

e{q) =2h- Jcosq. (4.3) 

Note that 1/(1 + 0^^"^^/-^) is nothing but the Fermi distribution function. Compared with the zero 
temperature case, we can see that the integration range in each matrix element is replaced as 

/ dq^ p{q)dq. (4.4) 

Jkp Jo 

However, the asymptotic behavior of eq. (4.2) is quite different from eq. (2.10) . In fact, this time 
we can apply Szego's theorem to the Toeplitz determinant like (4.2) (see Chapter 10 in ref. 31). 
The theorem tells us 

oo 

\nP{n) = npo + '^kpkP-k + o{l), (4.5) 



k=l 



where 



lnp{q) = £ Pfce'^-?. (4.6) 



In particular, we have 



fe=— oo 



1 /"^'^ 

^"^^'2^1 ln(l + e--('?)/^)dg, (4.7) 



Now recall that the free energy per site for the isotropic XY model is given hj^^'^^^ 

fiT,h) = -h--j^ ln(l + e-^(<')/^)d9. 
Then we can conclude that P{n) decays exponentially at finite temperature as 



T 

2'', r fJcosq-2h 
m ' — 



(4.8) 



P(n) ~ c(T, h) exp y-jj- (4-9) 

where the inverse of the correlation length and the prefactor c{T, h) are given by 

-1 f{T,h) + h 



dg, (4.10) 



1 + exp 

c(r, h) = exp ( ^ fcpfcP-fe) , (4.11) 
respectively. Note a formula 
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which can be found on page 112 of ref. 32. It is useful for numerical evaluation of c(T, h) (4.11). 

At finite temperature, it is known that the correlation functions such as {S^Sq) and {{S^Sq)) 
also decay exponentially. The correlation lengths of these correlation functions have been calculated 
analytically in refs. 15 and 33-35 as 



1 



2tt 



J cos q — 2h 



2T 



= smh ( I + smh 



(4.13) 
(4.14) 



In Fig. 3 and Fig. 4, we show some plots of -P(n) calculated from eq. (4.2) for T = 0.05 and 
r = 1. The dotted lines are the asymptotic functions (4.9). 



10 
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Fig. 3. P{n) at a low temperature : J = 1, T = 0.05. 

We observe the asymptotic behaviors of P{n) arc actually described by eq. (4.9). Particularly 

in the high temperature case (T = 1), even P[n) for small n are well fitted by eq. (4.9). At low 

temperature (T = 0.05), P[n) deviates from eq. (4.9) for small n. At this region, P(n) may behave 

more like Gaussian as in the zero temperature case. 

It is known that, when \h\ < h^, the low temperature expansion of the free-energy (4.8) is given 
byn,36) 



fiT,h) = e-^ + 0{T% 



(4.15) 



where e is the ground state energy (2.8) and Vg is the velocity of a low energy excitation 

Vs = Jcos-^ {2h/J) . (4.16) 
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Fig. 4. P{n) at a high temperature : J = 1, T = 1. 



Prom eq. (4.15), we can obtain the low temperature expansion of ^ ^ as 

1 

+ 



(Jsin [cos"^ C^h/J)] - Ihcos'^ i'^h/J)} 
ttT 



+ 0(T^). 



(4.17) 



6Jcos-^{2h/J) 

Thus we find the inverse of the correlation length diverges when T ^ as long as \h\ < he- 
That is, the correlation length ^ approaches to zero in the low temperature limit, which signals 
a crossover of the exponential decay to a stronger (Gaussian) decay in §2 and §3. This is in a 
remarkable contrast to the low temperature behaviors of ^xx (4-13) and ^zz (4-14). They behave as 
^ 1/T in the low temperature limit, which results from a crossover to a weaker (algebraic) decay. 

When the magnetic field is critical, i.e., h = ±/ic {he = ^/2), the low temperature expansion 
of the free-energy takes a differnt form,^^^ 



fiT,h) = -h,- 



1 



1 



J_\ ^3/2 

V2j 



+ o(t2). 



(4.18) 



Accordingly, the low temperature behavior of ^ ^ at h = ±hc is given by 



V2^ V2 



i- — )Vt + o{t) 



2hc 
T 



1 



V2TrJ 



c 



(for h = he 
+ 0{T), (for h = - 



(4.19) 



he) 



respectively. Thus we find particularly that goes to zero as ~ \/T when h = he- 

In Fig. 5, we plot the temperature dependences of We observe the low temperature 
behaviors of are actually described by eqs. (4.17) and (4.19). We also plot the prefactor c{T, h) 
in Fig. 6. When h < he, the prefactor c(T, h) diverges in the low temperature limit. 



Masahiro Shiroishi, Minoru Takahashi and Yoshihiro Nishiyama 




Fig. 6. Temperature dependence of the prefactor c{T,h) : J = 1. 
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In the high temperature hmit (T — ^ oo), we have 

P(n)^Qy, (4.20) 

since 

r'^ln(2), c(r,^)^l. (4.21) 

Alternatively, we can regard the magnetic field h as a function of the temperature and fix the 
magnetization 

In this case, since 
we have 

P{n)^(^^ + m^ =P{ir, (4.24) 

as r — oo. 

§5. Numerical Simulation by QMC Method 

In this section, we perform numerical simulations in order to confirm the analytical theory for 
P(n) at finite temperatures described in the previous section. Here we employ the quantum Monte- 
Carlo (QMC) method.^^) For finite-temperature calculations, the QMC method is particularly of 
use. In fact, recently, there have been proposed a number of substantial improvements: In ref. 38, 
an algorithm involving infinite Trotter-decomposition number is presented; that is, the continuous- 
time algorithm. Employing this algorithm, the authors had demonstrated that the simulation result 
is completely free from the Trotter-decomposition error. In addition, with the use of global update 
algorithm postulated in refs. 39-41, the auto-correlation time is reduced to considerable extent so 
that we can avoid wasting Monte-Carlo steps (the so-called critical slowing down). Hence, Monte- 
Carlo simulation combined with these techniques would be promising for studying long-range form 
of P{n) precisely. 

We treated systems with size N = 128, and imposed the periodic boundary condition. We 
performed five-million Monte-Carlo steps initiated by 0.5 million steps for reaching thermal equi- 
librium. P(n) is measured and averaged over the five-million Monte-Carlo steps. The results for 
T = 1 are plotted in Fig. 7. We see that the data are governed by the analytical asymptotic function 
postulated in §4 (shown by the dotted lines). We have performed the simulations for several other 
temperatures and magnetic fields. The obtained data exhibit good coincidence with the analytical 
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ones unless the temperature is exceedingly low (T < 0.1). In this way we could confirm the validity 
of the asymptotic formula (4.9) numerically at least for moderate temperatures. Unfortunately 
when the temperature is very low, the QMC simulation does not work very efficiently so that we 
could not reproduce the data such as at T = 0.05 yet. 
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Fig. 7. P(n) calculated by QMC method : J = 1,T = 1. 



§6. Summary and Discussion 

In this paper we have studied P{n) for the isotropic XY model in detail. Especially we have 
obtained the expression of P{n) in terms of Toeplitz determinant. On the basis of the expression, 
we could calculate the exact values of P(n) for small n and also get the asymptotic expression as 
n — oo. 

At zero temperature, Widom's theorem for Toeplitz determinant gives us the complete asymp- 
totic function directly. It agrees and completes the previous known results on the asymptotics of 
P(n). The asymptotic form of P{n) is given by 

P(n) ~ a-"' n", (6.1) 

with 

a = 1.4142, a = -0.25, (6.2) 

for zero magnetic field. For small n, we have also calculated P{n) by means of the DMRG method. 
In general, DMRG is not so efficient to investigate system at criticality such as eq. (2.1), because 
any correlation evaluated with the method, in principle, decays exponentially at large distances. 
However, we have found that as for P{n), which decays like a Gaussian, the method yields precise 
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result as well as correct asymptotic form. We could also estimate the speed of Gaussian decay a, 
etc., numerically. The method will be applied to other models in the forthcoming papers. 

We could obtain P(n) in terms of Toeplitz determinant also at finite temperature. This time 
from Szego's theorem, we have shown analytically that P{n) decays exponentially as n — oo, 

P(„)~exp(I±^>i. (6.3) 

We have confirmed this formula independently by means of the QMC method. 

It is Boos and Korepin^^ who first obtained the asymptotic formula (6.3) for any A when h = 0. 
They derived the formula from an observation of the definition of the thermodynamics, which can 
be generalized to the finite h case straightforwardly.^^) In the forthcoming papers, we will plan to 
confirm (6.3) for other models by use of the QMC method. It may be also possible to discuss P{n) 
at finite temperature from the point of view of the quantum transfer matrix method. ^^"^^'^^"^^^ 
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Appendix A: Equivalence of eq. (2.13) and eq. (2.10) 

Here we will show (2.15), i.e., 

-J ^det VeK-^'-^-)^- 

n\Jk, 27r Jk, 27r J,^^^/ 



is transformed to the Toeplitz determinant (2.10). Actually from 



det 



n— 1 

gi(lj!-l?m)j 

3=0 



det [e'«'^] det 



l,m=l 



= det 



.1=1 



n-1 



j,k=0 



we can see (A-1) is equal to 



JkF 2'"" Jkp 



det 



n— 1 



j,k=0 



det 



27r — fen 



27r 



n-1 



J j,k=0 



The last expression is identical to the Toeplitz determinant (2.10). 



(A-1) 



(A-2) 



(A-3) 
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Appendix B: Equivalence of eq. (2.17) and eq. (2.10) 

We will show the multiple integral formula (2.17) is reduced to the Toeplitz determinant (2.10). 
Here we assume kp > — for simplicity. 



Prom the relation 



TT 

sinh(x — — i) = — icoshx, (B-1) 



the integral formula (2.17) is written as 



, , . n(n-l) n(n+l) dAl dA„ 



-TT . , -A- siny-^(A,- - 5i)sinh"-^(A,- + fi) , , 

n 2 sinh(Aa - A,) J] j',,. (B-2) 



^ ■ ■ " cosh"2A,- 

a>b j=l ■' 



Furthermore by using the relations 



]^ 2sinh(A„ - Afo) = J| (^e^-"^" - e^f-"^" 

a>b a>h 



= det fe(2^-"-i)^^T , (B-3) 
L J j,k=i 

and 

cosh2Aj = 2sinh ^A^ - ^i^ sinh (^Xj + ^ij , (B-4) 
we can find that eq. (B-2) is represented by the determinant, 

n(n—l) 7i(n+l) 

P(n) = (-i) — 2— det[Ij,]l,^, (B-5) 

where 

dA 2e2^^ 

~ 2^(e2^-i)"-J-i(e2^ + i)^ ' ^ ^ ^ 

Now we introduce a change of integration variable from A to x defined by 

e2A = 1. (B.7) 

e^^ + i ^ ^ 

Then after some simple calculations, one can find the determinant representation (B-2) becomes 

n(l — n) ~ 

P(n) = (-2i)— det[I,fe]^,,=i, (B-8) 



where 

f27r— fep 
'fcp 



/■2-K-kF J 

ij^ = J ^(e^^ + (ei^_i)fe-V(J-")^. (B-9) 
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Recall that the relation between A and fcp is given in eq. (2.19), 

„-2A • ^ ^ 
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(B-10) 



Below we will show eq. (B-8) is identical to eq. (2.10). First we subtract the (n — l)-th column 
from the n-th column. Then the integrand of the n-th column can be replaced as 



(e'^ - i)"-ieHJ-»)^ ^ -2z(e'^ - i)"-2eiO-n)a;_ 



(B-ll) 



Similarly, wc can subtract the (A; — l)-tli column from the k-th column subsequently until k = 2. 
This procedure allows us to replace of the matrix elements for A; > 2 as , 



.2i / — (e^^ + i)"-'=(e^^ - i)fc-2ei0-")^. 
JkTP 27r 



(B-12) 



The above transformation yields a factor (— 2i)"~^, which cancels parts of the prefactor (B-8). 

Next we apply the same transformation to the n-th column again and subsequently to the A;-th 
column down to A; = 3. This procedure reduces the order of the factor (e^^ — i) in the integrand by 
1 and at the same time generates a factor (— 2i)"~2 Repeating the similar procedure, we finally 
arrive at 



P(n) = det 



27^-% Ax . , ■/ ■ N 



(B-13) 



Expanding (e^^ + i)" with respect to e'^, we can further simplify eq. (B-13) as 



P(n) = det 
= det 



Ap 27r 

/ 

Jk 



{n-k)x^i{j-n)x 



j,k=l 



r2iT-kF 



dx 



\{j-k)x 



(B-14) 



3,k=l 



The last expression is nothing but the Toeplitz determinant (2.10). 



1) V. E. Korepin, A. G. Izergin and N.M. Bogoliubov: Quantum Inverse Scattering Method and Correlation Functions 
(Cambridge University Press, Cambridge, 1993). 

2) M. Jimbo and T. Miwa: Algebraic analysis of solvable lattice models (AMS, 1995). 

3) V. E. Korepin, A. G. Izergin, F. H. L. Essler, D. B. Uglov: Phys. Lett. A 190 (1994) 182. 

4) F. H. L. Essler, H. Frahm, A. G. Izergin and V. E. Korepin: Commun. Math. Phys. 174 (1995) 191. 

5) F. H. L. Essler, H. Frahm, A. R. Its, V. E. Korepin: Nucl Phys B 446 (1995) 448. 

6) N. Kitanine, J. M. Maillet, V. Terras: Nucl. Phys. B 567 (2000) 554. 

7) A. V. Razmov and Yu. G. Stroganov: cond-mat/0012141, cond-mat/0102247. 

8) H. E. Boos and V. E. Korepin: hep-th/0104088, hep-th/0105144. 

9) M. Takahashi: Thermodynamics of One- Dimensional Solvable Models (Cambridge University Press, Cambridge, 
1999). 

10) E. H. Lieb, T. Schultz and D. Mattis: A. Phys. (N.Y.) 16 (1961) 417. 

11) S. Katsura: Phys. Rev. 127 (1962) 1508. 



20 Masahiro Shiroishi, Minoru Takahashi and Yoshihiro Nishiyama 

12) T. Niemeijer: Physica 36 (1967) 377. 

13) B. M. McCoy: Phys. Rev. 173 (1968) 531. 

14) S. Katsura, T. Horiguchi and M. Suzuki: Physica 46 (1970) 67. 

15) E. Barouch and B. M. McCoy: Phys. Rev. A 3 (1971) 786, 2137. 

16) H. G. Vaidya and C. A. Tracy: Physica. 92A (1978) 1. 

17) T. Tonegawa: SoUd State Commun. 40 (1981) 983. 

18) F. Colomo, A. G. Izergin, V. E. Korepin and V. Tognetti: Phys. Lett. A 169 (1992) 243. 

19) A. R. Its, A. G. Izergin, V. E. Korepin and N. A. Slavnov: Phys. Rev. Lett. 70 (1993) 1704. 

20) J. Stolze, A. Noppert and G. MiiUer: Phys. Rev. B 52 (1995) 4319. 

21) M. L. Mchta: Random Matrices, Second Edition, (Academic Press, 1991). 

22) C. A. Tracy and H. Widom: Comm. Math. Phys. 163 (1994) 33. 

23) P. A. Deift, A. R. Its and X. Zhou: Ann. Math. 146 (1997) 149. 

24) P. Deift: Orthogonal Polynomials and. Random Matrices: Riem,ann-Hilbert Approach (AMS, 2000). 

25) M. Jimbo and T. Miwa: J. Phys. A 29 (1996) 2923. 

26) Y. Fujii: J. Phys. Soc. Jpn. 69 (2000) 3143. 

27) H. Widom: Indiana Univ. Math. J. 21 (1971) 277. 

28) S. R. White: Phys. Rev. Lett. 69, (1992) 2863. 

29) S. R. White: Phys. Rev. B 48, (1993) 10345. 

30) Density-Matrix Renormalization — A New Numerical Method in Physics, ed. I. Peschel, X. Wang, M. Kaulke 
and K. Hallberg (Springer- Verlag, Berlin, Heidelberg, 1999). 

31) T. T. Wu and B. M. McCoy: The two-dim.ensional Ising m,odel (Harvard University Press, Cambride, 1973). 

32) M. L. Mehta: Matrix Theory (Editions de Physique, Orsay, France, 1989). 

33) M. Inoue and M. Suzuki: Prog. Theor. Phys. 79 (1988) 645. 

34) M. Takahashi: Phys. Rev. B 44, (1991) 12382. 

35) A. Kuniba, A. Sakai and J. Suzuki: Nucl. Phys. B 525, (1998) 597. 

36) M. Takahashi: Prog. Theor. Phys. 50 (1973) 1519. 

37) M. Suzuki: Prog. Theor. Phys. 56 (1976) 1454. 

38) B. B. Beard and U.-J. Wiese: Phys. Rev. Lett. 77 (1996) 5130. 

39) H. G. Evertz, G. Lana and M. Marcu: Phys. Rev. Lett. 70 (1993) 875. 

40) U.-J. Wiese and H.-P. Ying: Z. Phys. B 93 (1994) 147. 

41) N. Kawashima and J. E. Gubernatis: Phys. Rev. Lett. 73, (1994) 1295. 

42) S. Ostlund and S. Rommer: Phys. Rev. Lett. 75 (1995) 3537. 

43) V. E. Korepin: private communication. 

44) M. Suzuki and M. Inoue: Prog. Theor. Phys. 78 (1987) 787. 

45) T. Koma: Prog. Theor. Phys. 78 (1987) 1213. 

46) T. Koma: Prog. Theor. Phys. 81 (1989) 783. 

47) J. Suzuki, Y. Akutsu and M. Wadati: J. Phys. Soc. Jpn. 59 (1990) 2667. 

48) M. Takahashi: Phys. Rev. B 43 (1991) 5788. 

49) A. Kliimper: Z. Phys. B 91 (1993) 507. 



